On the size and shape of excluded volume polymers confined between parallel plates 
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A number of recent experiments fiave provided detailed observations of the configurations of long 
DNA strands under nano-to- micrometer sized confinement. We therefore revisit the problem of an 
excluded volume polymer chain confined between two parallel plates with varying plate separation. 
We show that the non-monotonic behavior of the overall size of the chain as a function of plate- 
separation, seen in computer simulations and reproduced by earlier theories, can already be predicted 
on the basis of scaling arguments. However, the behavior of the size in a plane parallel to the 
plates, a quantity observed in recent experiments, is predicted to be monotonic, in contrast to the 
experimental findings. We analyze this problem in depth with a mean-field approach that maps the 
confined polymer onto an anisotropic Gaussian chain, which allows the size of the polymer to be 
determined separately in the confined and unconfined directions. The theory allows the analytical 
construction of a smooth cross-over between the small plate-separation de Gennes regime and the 
large plate-separation Flory regime. The results show good agreement with Langevin dynamics 
simulations, and confirm the scaling predictions. 

PACS numbers; 82.35.Lr,36.20.Ey, 87.15.-v 
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I. INTRODUCTION 



With the typical length scales in cells ranging from 
tens of nanometers to tens of microns, their polymeric 
constituents are often spatially constrained. Some of the 
prominent examples are the cytoskeletal filamentous pro- 
tein aggregates microtubules and F-actin, which can have 
lengths up to tens of microns. Recent studies on micro- 
tubules in plant cells DNA packaging in viral cap- 
sids 2], DNA segregation in bacterial cell division 
focused on properties of biopolymers strongly influenced 
by confining geometries. This has naturally led to in- 
creased interest in understanding the physics of strongly 
confined polymers. Among all the polymeric constituents 
of cell, DNA with its relatively low bending stiffness in 
comparison to its extremely long length has a special 
place as it can under most circumstances be described as 
an classical excluded volume polymer Q- Important ex- 
amples of confined DNA are (i) chromosomal DNA that 
can have bare lengths up to centimeters trapped inside 
cell nucleus of size in the range of 1 - 10/im, (ii) bacte- 
rial DNA of lengths between 0.1 - 100/xm trapped within 
a small bacterial volume which for the case of E.coli is 
about 2/zm long and 0.5/zm in diameter, (iii) mitochon- 
drial DNA of length about 5/im for human beings with 
the available mitochondrial size between 0.5 - 10/im. The 
combined effects of confinement, self-avoidance and en- 
tropic forces act together in deciding the structure and 
properties of confined polymers. 

A number of single molecule fluorescence microscopy 
studies have recently investigated the structure and 
dynamics of confined DNA trapped between parallel 



'Electronic address; chaudhuri @amolf.nll 
t Electronic address; mulder@amolf.nll 



plates [§-[Tl| or within cylindrical geometries [l2l - [l3 | with 
the confining dimensions typically of the order of, or 
smaller than, the bulk radius of gyration of the polymer. 
In particular Bonthuis et al. [10] have recently observed 
DNA confined to parallel-plate nanochannels, showing 
among others an interesting non-monotonic variation of 
the two-dimensional (2d) projected radius of gyration 
with the distance between channel walls. 

Theoretical analyses of confined excluded volume 
chains already have a long history, encompassing a num- 
ber of different techniques and approaches such as scaling 
theories [15,, JJJ , renormalization group methods [17,, J^] , 
computer simulations @, IT9l - l22| . and mean field the- 
ory [23| . The qualitative picture that emerges from these 
analyses is simple: With increased degree of confinement, 
the polymer size shrinks in the confining direction(s), and 
expands in the non-confined direction(s). The expansion 
in the non-confined directions is due to the excluded vol- 
ume repulsion between different polymer segments. Ac- 
cording to de Gennes' scaling theory [l^, this expansion 
in the strongly confined regime has a power-law depen- 
dence on the length-scale of the confinement. On the 
other hand, in the limit of very weak or no confinement 
the polymer is expected to behave as a free excluded vol- 
ume polymer approximately obeying the Flory relation- 
ship between polymer length and spatial size in three 
dimensions (3d) Q- 

Monte-Carlo (MC) simulations of self-avoiding lattice 
random walks trapped within refiecting walls have vali- 
dated the de Gennes scaling predictions [l^, [l^ . These 
simulations also showed a non-monotonic variation of the 
mean size of the polymer, measured as the full 3d radius 
of gyration, as a function of the inter- wall separation . 
This non-monotonicity in the radius of gyration was later 
also cap tured by a mean field theory (MFT) calcula- 
tion [231 . However, the experiments reported in Ref. [To[ 
found non-monotonicity in the projected 2d size of the 



polymer, rather than its 3d average. To correctly inter- 
pret these results we need a theoretical approach that ex- 
plicitly allows us to consider the behavior of the polymer 
in the confined and non-confined directions separately. 

Here we revisit the problem of an excluded volume 
polymer confined between two parallel plates with the 
latter desideratum in mind. By reviewing the extant 
scaling approaches, we show that the non-monotonicity 
of the 3d radius of gyration with plate-separation already 
follows from a scaling analysis if due care is given to the 
details of its application. We then construct a mean-field 
theory (MFT) for the confined polymer that goes be- 
yond the existing theory, in that it explicitly takes into 
account the symmetry breaking induced by the confine- 
ment in the construction of the non-interacting reference 
polymer. This approach allows the analytical determi- 
nation of the prefactors of the scaling relations in the 
limits of strong and weak confinement respectively. The 
results of this approach are then validated by an explicit 
off-lattice simulation of a confined polymer. Our main 
result is that the size of the polymer transverse to the 
confining direction, i.e., parallel to the confining planes, 
grows monotonically with decreasing separation between 
the plates. The non-monotonicity of the 3d radius of gy- 
ration therefore is due to the fact that, with decreasing 
plate-separation, initially the polymer-size in the confin- 
ing direction decreases more rapidly than that it grows 
in the unconfined directions, and eventually the latter 
becomes dominant. 

The paper consists of three main sections. We be- 
gin with by describing the existing scaling approaches in 
Sec. ini In Sec. IIIII we present our variational MFT cal- 
culation and its predictions. In Sec. IIVI we describe the 
Langevin dynamics simulation scheme, and compare the 
results with our MFT predictions. Finally we conclude in 
Sec.|V]by summarizing our results and providing an out- 
look. A number of appendices collect some of the more 
technical material. 



II. SCALING THEORIES 

In this section, we briefly review the existing scaling 
theories for a polymer confined between parallel plates. 
We show that the non-monotonicity of the mean polymer 
size as a function of the plate separation, as observed in 
earlier simulations [l^, is in fact intrinsic to de Gennes 
scaling predictions. 

The first scaling theory of self-avoiding polymer be- 
havior is due to Flory. Because of their simplicity Flory 
like arguments have been extended to other cases, e.g., 
semifiexible polymers |2J]. Within Flory scaling the- 
ory the free energy of an excluded volume polymer in 
d-dimensions can be written as '23 , 
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FIG. 1: (Color online) Schematic diagram of the non- 
monotonic behavior of average polymer size Rt^t ~ aiD^ + 
a2D^^^^ (dark thick line). In the limit of small D {< L con- 
tour length), the combination of and D^^^^ behavior of 
the polymer size in the confining and orthogonal directions, 
respectively, leads to non-monotonicity with a clear minimum 
in the average polymer size Rt^^ . At large D (> L), R^^^ 
crosses over to the Flory regime and becomes independent of 
D. 



structed of N segments, each of size I, (3 = l/fcsT with 
ks the Boltzmann constant and T the ambient temper- 
ature. The first term describes the entropic elasticity of 
a free chain, and the second term describes the repulsive 
interaction between different segments. Minimizing the 
free energy with respect to R one obtains the Flory esti- 
mate of the equilibrium polymer size R ~ N^/('^+'^) I. As 
de Gennes pointed out [7,], the success of Flory's theory 
relies on remarkable cancellation of overestimates in both 
the terms. For polymers confined between parallel plates, 
a Flory like argument predicts an equilibrium polymer- 
size in the unconfined directions that agrees with de 
Gennes scaling (which we discuss next), but the free en- 
ergy it predicts is not extensive in N |22| . 

The de Gennes scaling arguments for the size of a poly- 
mer confined between parallel plates can be described 
as follows [l^ . If the plate-separation D < L — N I 
the polymer contour length, the size of the polymer in 
the confining direction (perpendicular to the confining 
plates), i?_L, becomes limited by the confinement, so that 
R^ ~ D. The size of the polymer in the unconfined di- 
rections (parallel to the confining plates), should be 
entirely decided by the following two length scales: The 
3d Flory size of a polymer R ^ N^^^l, and the plate- 
separation D. Thus = R(f){R/D), so that (/)(x) = 1 
when a; — )■ (3d Flory regime) and (/)(a;) ~ x'^ when cc > 1 
(de Gennes regime). In the 2d-limit (x ^ 1) Flory scaling 
requires R\\ ~ N'^^^l. This ensures q = 1/4 and there- 
fore ~ N^/^{l/Dfl^l. This power-law divergence of 
with decreasing D was observed in MC simulations 

19 , 2C I , and also in 
IflEI. The excess 
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where R is the end-to-end distance of the polymer con- free energy measured from the state of an unconfined self- 
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avoiding polymer can be expressed [T^ as /3F = ip{R/D), 
so that ^(x) = for X — >■ and i^{x) = in the other 
limit of small D. Demanding that the excess free energy 
has to be an extensive function of N one obtains p — 5/3 
and PF - N{l/Df/^. 

We notice that, in terms of plate-separation D de 
Gennes scaling predicts two completely opposite behav- 
iors for Ra^^ D and i?,, - N^/^{l/Dy/^ I. While Rjl de- 
creases, _R|| increases with decreasing D. Therefore, the 
mean size of the polymer, averaged over all the three di- 
rections, — R\ + 2i?| ~ aiD'^ -\-a2D~^/'^ (oi, a2 are 
positive constants) is expected to vary non-monotonically 
with D. The R^^^ ~ D^^^^ behavior at very small D 
crosses over to R^^^ ^ at moderately large D via 
a minimum in i^joj. This fact is illustrated in Fig. [T] 
and was observed in previous simulations [loj . At large 
D L), the scaling function (j){R/D) = 1 and therefore 
the average size R^^f. becomes independent of D (Fig. [ij . 
It should be noted here that this behavior is not the same 
as the non-monotonicity observed in experimental mea- 
surement of the parallel component of the radius of gyra- 
tion of a DNA confined between parallel plates 10, 11]. 
de Gennes scaling predicts a clear monotonic increase of 
the component i?|| with shrinking D. 

It is instructive to compare this behavior with that of 
an ideal chain confined between parallel plates. For an 
ideal chain, in absence of any interaction, the impact of 
the confinement remains restricted to the confining direc- 
tion only. Thus R± ^ D for small D and R± ~ N^^^l in 
the bulk limit. However, _R|| ~ N^^^l is obeyed indepen- 
dent ofD. Therefore the overall size is R^^^ ~ a3D^+2_R^ 
with 03 a positive constant. Starting from the bulk value, 
with shrinking D, R^^^. reduces to ultimately saturate at 
2i?^. The initial reduction of overall size is thus common 
to both ideal polymers and excluded volume polymers. 
However, the expansion of the overall polymer size with 
further reduction of D for excluded volume polymers is 
a signature of the repulsion between polymer segments. 



III. MEAN FIELD THEORY 

A. Overview 

The starting point for our MFT is a description of the 
polymer as a space curve r(s), in which s is the arc length 
parameter with domain < s < L. The energetics of the 
excluded volume (i.e., self-avoiding) polymer is given by 
the Edwards Hamiltonian [2BI 

+ ^bj^ dsj^ d.'<5(r(s)-r(s')). (1) 

The parameter b sets the strength of the repulsive inter- 
segment interactions. 



Edwards and Singh [26*1 first proposed to determine the 
size of a free self-avoiding polymer by mapping it onto an 
isotropic Gaussian chain with reference Hamiltonian 

and then choosing the effective segment size Iq self- 
consistently by requiring the first order variation of the 
mean-squared end-to-end distance {R"^) in the interaction 
parameter b to vanish. This approach relatively simply 
reproduces the Flory result R ~ N's for the scaling of 
polymer size with polymer length. This same idea was 
later applied by Thirumalai and co-workers [2^, [13] to 
confined polymers by first imposing the confinement on 
the isotropic Gaussian reference chain and then perform- 
ing the self-consistent calculation of the effective segment 
length. 

One can, however, pose the question whether an 
isotropic ideal reference chain is appropriate to the sit- 
uation where the presence of the boundaries already ex- 
plicitly breaks the bulk rotational symmetry. Here we 
propose to map the confined polymer onto an anisotropic 
reference chain, composed of bonds whose length depends 
on their absolute orientation with respect to the sym- 
metry axes of the confining geometry. This anisotropic 
Gaussian reference- chain is described by the Hamiltonian 

= m\\ +pn±. (2) 

The effective segment lengths ?|| and l± are now the pa- 
rameters that have to be determined self-consistently. 
Apart from the fact that this choice is physically plausible 
given the geometry of the system, it has two other po- 
tential advantages. It more explicitly yields independent 
predictions for the behavior of the polymer in the paral- 
lel and transverse directions with respect to the confine- 
ment, and generically multi-parameter mean field theo- 
ries can be expected to yield better approximations to 
the underlying physics than single parameter ones. 

We assume the z-axis of our coordinate system to be 
perpendicular to the confining planes. For convenience, 
we consider the lateral extent of the system to be finite, 
defined by a length W ^ D, but freely take the limit 
of an infinite extent wherever appropriate. The hard- 
wall confinement is implemented through the Dirichlet 
boundary condition by which the Green's function cor- 
responding to the reference Hamiltonian vanishes at the 
two confining walls. The full Hamiltonian can be ex- 
pressed as 

I3H = mo + /3A-H|i + l3An± + fSAHb (3) 



3 



where 

PAnb = l-b [ ds [ ds'(5(r(s) -r(s')) (6) 
2 Jo Jo 



where the normaUzation constant is given by the parti- 
tion function 

rfr(O) / dv{L) / I?[r(s)]exp{-/3H[r(s)]}. 

J Jr{0) 

(8) 

We now consider the components of the end-to-end 
separation R|| = r|| (L) — r|| (0) and R± = r± [L) — r± (0) 
and the averages and {R']_)- Expanding the Hamil- 
tonian Ji in the above expression up to hnear order in 
AH around the reference Hamiltonian T-Lq, we obtain 

(i?}j) = (i?f)^-<5|l+0(AH2), (9) 
{Rl)={Rl)^-S^+OiAn') (10) 

where 

<5|l = [{RpAH\i)o - (i?|)o(/3AH||)o} (11) 
+ {{RpAni_)o - (i?f)o(/3AHi)o} 
+ [{RpAnb)o - (i?|)o(/3AHb)o} 

S± = {(i?i/3AH||)o - (i?i)o(/3AH||)o} (12) 
+ {(i?l/3AHi)o - (i?i)o(/3AH±)o} 
+ {(i?i/3AH,)o - {Rl)oWAHb)o} 
= s{ + 6i+S'i. 

{. . .)o indicates a statistical average taken with respect 
to the reference Hamiltonian Hq. In the above expres- 
sions, (5||-, Sj^ and are the corrections due to the 
chain anisotropy, and (5|j and S'^ are the corrections from 
excluded volume interactions. The MFT approximation 
requires a choice of Z|| and l± so that {R'p = (i?|)o and 

(i?^) — {R\)o- Thus, solving = (5x = 0, one can ob- 
tain the effective segment lengths Z|| and l±_ in terms of 



are treated as perturbative corrections. 



Statistical averages with respect to the Hamiltonian H 
are defined through a path integral 



(7) 

I 

plate-separation D, contour length L and inter-segment 
interaction strength b. 

Ref. 1261] obtained Flory scaling within corrections up 
to linear order in PAI-L. The inclusion of higher order 
terms was shown to change the coefficients of scaling 
form, however, keeping them stable [26^. We expect the 
same to hold in the present problem as well [23], as the 
nature of interactions used here is exactly the same. 

We now briefly summarize the main results we ob- 
tained using this MFT. In the limit of very large plate- 
separations {D — > 00) the system becomes isotropic with 

1\\=1^ = Ib^ {blf'^'L^''' (13) 

and the mean squared end-to-end separation follows the 
relation 

= LIb - {blf'^L^'^ (14) 

which is independent of D and obeys Flory scaling. 
In the limit of strong confinement {D — >■ 0), 

/ / I 4&zJ\ , 4V2^ [bD^/^, , 

Using these, we obtain the size of the polymer in the 
confining direction 

with K = {\ / Q)l / D)"^ , and in the unconfined direc- 
tions 

^ ,/X(M)l/2/,3/27^-l/2. (18) 



{A[v{s)])^^ j dr(0) / dv{L) r''^^2?[r(s)]A[r(s)]exp{-/3H[r(.)]} 
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With shrinking plate-separation D the polymer gets com- 
pressed in the confining direction and expands in the free 
directions. Up to the leading order (in the small D limit) 
{RY) and ~ L^/^D-Vs obeying de Gennes 

scaling. Beyond the leading order, the expressions we 
found describe a smooth crossover to large D behavior 
(see Fig.E]). 

We now present a detailed derivation of the Eqs. [T^- 

m 



B. Formulation 

1. Green's function 

In order to determine all the above mentioned aver- 
ages, we require the Green's function corresponding to 
the anisotropic reference Hamiltonian Ho [Eq. ^] de- 
fined as 



where 



Z|| = J drii J drji G||(r||,i | rl|,0) = (25) 



and 



-KL{2n+l) 



= E 



(2n + 1)%2 

n— ^ 



(26) 



Remember that denotes the area of the unconfined 
xy-plane, and we freely take the limit W ^ oo wherever 
appropriate. 



3. End-to-end separation: Contribution from reference 
Hamiltonian 



nT=r(L) 

Go (r,i|r', 0) - / V[r{s)] exp {-/3Ho[r(s)]} . 

Jr'=r(0) 

(19) 

The Green's function obeys the following differential 
equation [28|], 



(20) 



Using the Dirichlet boundary condition in the confin- 
ing direction Go{z — 0, D) — and free boundary con- 
dition in the unconfined directions (along x and y axes), 
one obtains a variable separable form of the Green's func- 
tion (see Appendix-El 

Go-G,,(r,|,L|r;,,0) G^{r^, L \ r'^,0) (21) 



Using the above mentioned Green's function Go one 
can calculate the end-to-end separations 



/dr,|/rfrj| (r||-rj|)'G|| (r,, ,L|rj| , o) 



and 

{Ri)o 



3 " 



d r^ drl (r^ - 4)^ G^ (r^,£|rl, 0) 
/o^dr^/o^drlG^ (ri,i|rl,0) 



/dr||/dr|iG|| (r||,L|r|,,0) 



(27) 



2 ? 



(28) 



where, 



where. 



and 



G|i 



:(r||-r;,)' 



27rZ||L 



nTT 



D 



-Kn L 



(22) 



(23) 



1 - 



-K{2n+iyL 



% {2n+lY \ (271 1)'' 7r2 



^1 (2n)^ 



-K(2n)^L 



and 



(29) 



with K - (§) 



^. Partition function 



OO ^ 

^ = E — — 



,-K(2n+l)''L 



(30) 



The partition function corresponding to the refer- 
ence Hamiltonian "Hq can also be written in a variable- 
separable form 



Z = Jdr{L) y"dr(0) Go(r(L), L | r(0), 0) = Z,, Zj 



(24) 



4- End-to-end separation: Corrections from chain 
anisotropy 

We now calculate the first order contributions 
jjj , (5^, (5|[ and dj^ that are due to the anisotropy of the 
reference chain. Using the definition of the mean-squared 
end-to-end distance 
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/drii /dr|, (rii -r||)Vy;ff I?[r||(.)]exp{-/?H||[r(.)]} 
/ rfrii / rfrjj /;,"r;"i)^ (.)] exp [r(.)] } 



r 



along with the identity 



1 1 



/3^||[r(s)]=«|||(^y-^j /3AH||[r(s)], 
where = 3/2Z||, we obtain 



{(i?^/3AH||)o-(i?f)o(/3AH||)o} 



or, 5j[ 



1 1 



Lli 



Similarly, one can show that (with aj_ — ^1/21^^, 
3/1 l\ d 
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(^^)o = 



2 V ; /, 



2 \ ^ Z|| / da\\ 
^ 2\l L dai 



(31) 

(32) 
(33) 
(34) 



where the last derivative can be explicitly evaluated by 
considering Eq. (pSj) . This leads to the following expres- 
sion, 



The calculation of these two terms are lengthy but 
straightforward. Here, we only quote the final results 
deferring the details of the calculations to Appendix-lBl 



= -Z-i- 



M{nL,n',no;K) 

,n',no — 1 

X [2(5rio,nt + ^n',{no+nL)/2 ~ ^n' ,{no-nL) /2 

-S„>,{nL-n„)/2] X Jo(»^L,"o) (36) 



where. 



M{nL, n', no; n) = e""-^"^ x 

L nL 

ds' I ds 



(12 2 \ ' I 2 /2\ 



, (37) 



Join L, no) 



noULn^ 



:[l-(-l)"-][l-(-l)"»]; (38) 



and 



6 3 

11/2 __ K{nL,n',no;K) 



Z ' 2ttI 

S. 



■II 

' nL,n\no — l 
X i'^Sng^riL ^ ^n',{no+nL)/2 



(no-nL)/2 ~ On' , (nL-no)/2 



JinL,no) 



{R±)q t I \ 

^2 JQ\nL,no) 



(39) 



where, 



12^ 

OO 



1 1 

8 



(2n+ l)27r2 



-kL(2ji+1)^ 



(35) 



where and ^ are defined by Eq. ([^5]) and Eq. ([50)) re- 
spectively. 



5. End-to-end separation: Corrections from excluded 
volume interactions 



K(nL,n',no;K) = e"''^"^ lim / ds' / ds 

'^^oJo Js' 



(s - s') + a 



m 



1 



J(ni,no) = -^^[(-l)»-ni^2^2(l-(-l)"-)] 

X [l-(-l)""] ?_^(-l)«o(_i)«x. 



^[(-l)-.§.2 + 2(l-(_l)-)] 



X [l-(-l)"-]. 



(41) 



Now we evaluate and (5^^. Note that the pres- 
ence of (5(r(s) — i"(s')) in the self-avoidance interaction 
(HAHb [Eq. ^] couples the transverse and the longitu- 
dinal modes of the Green's function. This is the term 
through which becomes dependent on the plate- 

separation D. 



C. Results 

1. Large plate-separation limit: Z) — >■ oo 

Here we demonstrate that the MET scheme naturally 
leads to Flory scaling in the bulk limit of D — >■ cxo. This 
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limit is better appreciated, in a shifted coordinate sys- 
tem in which the confining hard wahs are at z = ±Z)/2. 
In this coordinate system, the transverse component of 
Green's function is. 



n, ev< 
+ n E 



nnz niTZ ,.„2r 

n, even 

D ^ cos cos— e—^ 
ji, odd 



In the continuum hmit of Z) — > cxi 

i{z-z'f 



exp 



2 Ia_L 



(42) 



(jjj is independent of D and is given by Eq.[nil The details 
of the calculation of , 5\ and in this bulk limit are 
presented in Appendix-ICl In this limit the equations, 

= 5] + (5|| = and 5^ = (5f + = imply, 



^3/2 



(43) 



and 



r/2 w^/^ 

7-3/2 



respectively. Eq. (gS]) and Eq. (gH) leads to Z|| = = ^s, 
the isotropy expected in the limit of Z? — > oo. This also 
implies 



Ib ^ {blf'L^"'. 



and 



(45) 
(46) 



= LIb - {blf'^^L^'^, 
i.e., Flory scaling in three dimensions. 

2. Narrow plate-separation limit: D — >■ 



In the limit of Z) — > 0, k = 1±{it/D)'^/6 — > oo and 
only the small n eigenfunctions contribute to the calcu- 
lation of perturbative corrections to end-to-end distance. 
The simplest approximation in this limit is the ground 
state approximation (only the minimum value of ti's con- 
tribute). 

Ground state approximation: We first calculate 
within the ground state approximation. Using uq = 
n' = UL = 1, we get M(1,1,1;k) = (L"^ /2) exp{-KL) 
and Z = W'^{8D/Tr'^)exp{-KL) which leads to (5|| = 

-{3b/An){L^/D) (using Eq. EH). Within MET, sat- 
isfies the condition Sh = + 6u = 0, which implies that, 



3 U hJ " 47r Z? ■ 



The solution of this equation gives 

which is Eq. 1161 Therefore, we find the relation given in 
Eq.im 




1 



-Ll 1 



95 L 
2^1 D 



In the last step, the diverging part of (Z?|)o in the limit 
of Z? — > is extracted. Thus, up to the leading order 
in the limit of Z? — )■ we find de Gennes scaling Rn = 



(Z?|) - {l/V2^){bl)^/^L^/*D-^^\ 

Now we calculate the perpendicular component (Z?^ 
within the ground state approximation = n' 
1. Thus Eg. [Ml reduces to 



no 



51 



96 



2nl\\ Z 



J(l,l) 



4 (Z?2 



D2 



In the expression of {R']_)o [Eq. ([28| ] we retain only the 
n — terms, within the ground state approximation, 
(44) to obtain {R\)o = 1 - S/tt^. This leads to J(l,l) - 
(4/7r2)(Z?i)o = and therefore, S'i = 0. It can be easily 
seen that within the ground state approximation (5^ = 
too. Thus /j_ remains indeterminate. We need to go up to 
the first excited state approximation (next higher values 
of n), to obtain the expression for ^j^. 

First excited state approimation: Within this approxi- 
mation, we can write Eq. ^\ as 



36 



27r/|i Z 



W^D^ [3K{1, 1, 1; k) Ti + 2K{1, 2, 1; k) Ti 



2K{2, 1, 2; k)T2 + 3K{2, 2, 2; k)T2] 



where. 



{R] 



L>2 



-•^0(1,1) 



27r2 



Z32 



'-M2,2) 



1 

2^' 





log 











Ti = J(l,l)- 
T2 = J(2,2)- 
We find. 



We observe that the leading order behavior of 
Zs'(1,1,1;k) ^ Lexp(— kL) apart from a weak log- 
arithmic divergence coming from the (5-function na- 
ture of the inter-segment repulsion (a — > 0). There- 
fore, Zir(l, 1, 1; k)Ti ^ icxp(— 4Ki). Similarly one 
can show that all the terms in the above expres- 
sion of 5']_ has the same leading order behavior. Us- 
ing Z = W"^ {8D/TT'^)exp{-KL) we find that, - 
(6LZ?/?||)exp(-3KL). 
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Within the first excited state approximation, we can 
write Eq. [35] as 

■""^ -3kL,2 t f ^ 1 



rj_ —OKL/lZ T 



Thus the mean field condition, 5± = 5j^ + S']_ — imphes, 











1 + 




~ 2 







-L — t: 



which is Eq. 1151 In the hmit of Z) ^ we can write 
/j^ ~ l-2bD/l\i ~ ;-(4n/2^/3)vV^-D3/^L-i/2. Thus the 
leading order behavior is ~ ^ with correction vanishing 
as £)3/2. Using this Z_l in k = {1/6)1±{tt/D)^, we find 



2 



1 



-3nL 



which is Eq. [TTl The above expression means that, the 
leading order behavior (at smallest D values) of 



{Rl) = 



1 - 



has corrections in the next order (at larger D), which has 
a complicated functional dependence on inter-segment 
interaction strength 6, plate-separation D and polymer 
contour- length L. 

Before ending this section, we briefly discuss the pure 
2d limit of the above-mentioned calculation. Variational 
MET calculation for pure 2d leads to the mean squared 
end-to-end distance 



(47) 



which obeys 2d Flory scaling. In the above expression 6o 
measures the strength of the inter-segment repulsion and 
is dimensionless, in contrast to the same parameter in 3d 
denoted by b which has the dimension of length. This 
indicates that I is the only intrinsic microscopic length 
scale in the system. Requiring b = bgl, we find that (^|) 
reaches the 2d limit of (R^d) when D ^ I (see Appendix- 



3. Non-monotonicity in overall polymer-size: Crossover 

Clearly, the mean squared end-to-end distance 
(-^tot) = -I- (i?^) should show a minimum as a 

function of changing plate-separation D. This comes 
about, because of the completely opposing behaviors of 
the two components - while the size in the confining di- 
rection {R']_) shrinks with decreasing D, the polymer-size 
in the unconfined directions (^|) expand. Keeping up to 

the leading order behavior in (i?^) [= (1 — 8/n^)D^ /2] 
one can easily find the cross-over plate-separation (using 
Eqs.[l2]and[18l) 



Z?, = (27r)-i/M 1-4 {biy/'L'/'' 



-2/5 



(48) 



where (i?tot) reaches the minimum. With enhanced de- 
gree of confinement (reduced D) the total size of the poly- 
mer (i?5_) fli'st shrinks and when D < Dc the size starts 
to grow. 

Dc follows the same 3d Flory-scaling form as the size of 
an excluded volume polymer freely floating in 3d. Thus 
the position of this minimum scaled by the bulk polymer 
size (end-to-end distance or radius of gyration) Dc/L^^^ 
should be independent of polymer contour-length L. The 
MC simulation results shown in Fig.l of Ref.|19| corrob- 
orate this fact. 



IV. SIMULATIONS 

In the previous two sections we demonstrated scaling 
arguments and a mean field theory to obtain estimates for 
the size and shape of confined excluded volume polymers. 
To validate our mean field picture, we present the results 
of a full numerical simulation for such a system. 



A. Method 

We perform molecular dynamics simulations of a self 
avoiding bead-spring chain trapped between two soft- 
walls in presence of a Langevin heat bath. All the beads 
interact via the fully repulsive part of Lennard-Jonnes 
potential V[f{r) = 4e [{a/rY^ - {a/rf + 1/4], with a 
cut-off distance set to Vc = so that V['j'{r > Vc) = 

0. e is the strength and a is the range of the interaction. 
The bond between two neighboring beads are modeled 
by a shifted harmonic potential Vsp{r) — {A/2){r — cr)^ 
with spring constant A — lOQe/a^ such that the equilib- 
rium bond- length, in absence of any other forces, is cr. r 
denotes the center-to-center distance between a pair of 
beads. We assume that the polymer is confined within 
two parallel walls placed at z = —D/2,D/2. The repul- 
sive potential due to the walls is assumed to be the inte- 
grated and shifted Lennard- Jones potential VwaiiiSz) = 
e [{cr/Szy° - {a/Sz^ + g] with g = (5/2-l)(2/5)5/3 and 

a cut-off distance Zc — (5/2)^/^(7, so that V^aH(|<^z| > 
Zc) — 0. Here Sz is the z-separation of a bead from any 
one of the walls (top or bottom), e and a set the energy 
and length scales, respectively. We integrate Langevin 
equations using the velocity- Verlet algorithmic^ with a 
time-step St = O.Olr where r — cr-^/m/e is the charac- 
teristic time scale. We set the mass of each bead m = 1 
and isotropic friction coefficient 7 = 1/t. The Langevin 
thermostat [33 is used to keep the system at a constant 
temperature T = l.Oe/fcs. 



B. Results 

The quantity we are interested in is the equilibrium 
size and shape of the polymer. In the model we simulate. 



8 



100 





/ 

/ 

/ 
















a 

a 


x(c) - 


□ 



1 10 100 



FIG. 2: (Color online) Components of mean squared end-to- 
end distance (i??) as a function of plate-separation D. All 
lengths are expressed in units of a. We have plotted the 
components (X^)' = {X^) - Ay, (Y^)' = (F^) - Ay, and 
(Z^)' = (Z'^) + Ax where \\\ and Ax are offset polymer sizes 
(see main text). With decreasing D, {X^) and (F^) increases 
in the same manner while {Tp') decreases. All the components 
obey de Gennes scaling, up to the leading order: (X^)' — 
(y^)' = ^t/-D^''^ and (Z^)' = "iD^ where p. = 306 and 7 = 
0.12. The fitted values of offsets are Ay = 12.8 ± 4.2 and 
Ax = 1.4. At larger _D, (Z^)' shifts from the scaling form 
(Eq[T7| asx(-D) = 7!?^ - ciD^ exp(-C2/L'^ -hcs/D^/^) where 
ci = 0.1, C2 = 970 and C3 = 2.06 ± 0.36. The fitting error in 
all the parameters is less than 5%, unless specified otherwise. 



the polymer is made of i ~ 1 . . . . , beads with a mean 
contour length L — (N — 1) in units of a. We follow 
the mean squared end-to-end distance in the confining 
^-direction (Z^) = ((z^r — zi)^) as well as in the uncon- 
fined xy-planc {X^) = {{xn~xi)^), {Y^} = {{yN-yi f}- 
The time evolution of these quantities allows us to iden- 
tify the equilibration of the system. In order to test the 
validity of our simulation scheme, in a separate simula- 
tion of the bulk-system (using periodic boundary con- 
ditions in all three directions) we measured {R'^^^y-) — 
{X^) + (F^) -I- (Z^) as a function of polymer contour 
length L = 8, 16, 32, 64, 128, 256 to obtain Flory scal- 
ing (i?buik) ^ ^^^^ (data not shown). 

For the simulations of confined system we used an iV = 
64 bead polymer. The simulations were equilibrated for 
lO^T before collecting data over further 3 x IO^t. We used 
periodic boundary conditions in x and y directions with 
the lateral extent of the simulation box in these directions 
W = 2L. In Figl2]we have plotted the three components 
of mean squared end-to-end distance {X'^)' = {X'^) — A||, 

and Aj^ are offset polymer sizes which comes about for 
reasons discussed in the following. 

The mean squared end-to-end distance in the confining 
direction {Z"^)' shrinks with increasing degree of confine- 
ment (decreasing D) as jD^ (FigH]). Apart from an ad- 
ditive offset Aj^, this is consistent with the leading order 
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FIG. 3: (Color online) Mean squared end-to-end distance 
(-Rtot) ~ (X^) + {Y'^) + {Z^) as a function of plate-separation 
D. The line is a plot of 2(^i/D^''^ + Ay) + (7!)^ - Ax) where 
we used the same values of n, 7, Ay and Ax as in Fig[2] 

behavior predicted by our MFT [Eq. (fTT]) ] and de Gennes 
scaling. Due to a finite non-zero range of repulsion Zc 
coming from the soft-wall confinement (not incorporated 
in the MFT as in the MFT we assumed hard wall confine- 
ment for simplicity), simulated {Z'^) gets suppressed by 
an extra amount A^ . This indicates that A^ should be a 
function of Zc which vanishes as Zc — (hard wall limit). 
The fitting procedure in Figl2]gives Aj^ = 1.4±0.06 which 
is numerically indistinguishable from z^ — 1.36. At larger 
D, we find a saturation of (Z^) that obeys the functional 
form (FigEl) 

X(D) = jD^ - ciD^ exp{-C2/D^ + c^/D^'^). 

Notice that this form of x(-D) is obtained from the ex- 
pression of {R\) obtained from Eqs. pT)) and (fT5)) . 

FigO clearly shows that with increasing degree of con- 
finement, the components of mean squared end-to-end 
separation in the unconfined directions expands {X"^)' — 
(Y^y _ y^/£)i/2^ Leaving out the offset value Ay, this 
obeys dc Gennes scaling. Note that an offset like Ay was 
expected from our MFT [see the first line of Eq. ^TE\i ]. 
Further, the excess compression of (Z^) by A_l requires 
an extra bit of expansion in the parallel components. Ay 
contains this contribution too. At large enough plate- 
separations (X^), (Y'^) saturate to their bulk value. 

The total mean squared end-to-end distance of the 
confined-polymer (i?^^^) = (X^) + (Y^) + (Z^) shows a 
non-monotonic dependence on D (Figl3|). (i?tot); start- 
ing from its bulk value at large D, first reduces to reach 
a minimum and then expands as we reduce the plate sep- 
aration. The initial decrease in size is mainly governed 
by the shrinkage in (Z^). At small D values, the large 
increase of the polymer-size in the unconfined directions 
(AT^) and (Y'^) takes over. The cross-over between these 
two different behaviors takes place at the minimum of 
(Rtot)^ a feature predicted by the MFT [Eq. (gS])] as well 
as by de Gennes scaling (SecHII- Fig [3] shows that the 
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scaling forms added with the offset parameters in the 
three components of the mean squared end-to-end dis- 
tance does capture the non-monotonicity as well as the 
approximate position of the minimum in (R'^^^} 



V. DISCUSSION AND OUTLOOK 

In this paper we have presented a mean field approach 
to calculate the mean squared end-to-end distance of 
a self-avoiding polymer confined within a parallel-plate 
geometry. The method allowed us to calculate all the 
components of this quantity separately as a function of 
plate-separation. Up to the leading order calculation, we 
recovered Flory scaling in the limit of large plate sepa- 
rations (bulk limit), and de Gennes scaling in the limit 
of small plate separations. The next to leading order 
correction allows a smooth transition from de Gennes 
regime towards Flory regime. We believe that higher or- 
der perturbation calculations using our MET framework 
will further bridge the gap between the two regimes. We 
noted that the non-monotonicity of the overall polymer 
size as a function of the increased degree of confinement, 
as clearly captured by our MET, was already inherent to 
de Gennes scaling. The numerical results from Langevin 
dynamics simulations showed good agreement with our 
MET predictions. 

It is, however, hard to make quantitative comparisons 
between the theoretical calculations and the simulations 
for quantities like actual prefactors of scaling laws. There 
are two reasons behind this. Eirst, these non-universal 
quantities depend on the detailed nature and strength of 
the interactions. Secondly, as was already discussed in 
Ref . [1^ , the effective values of the prefactors in the scal- 
ing forms depend on up to which order the perturbative 
calculations are performed in the MET. 

Before we end the discussion, we would like to em- 
phasize the simple fact that with decreasing plate- 
separations, the polymer size in the confining direction 
can only shrink, while it expands in the other uncon- 
fined directions. In a recent experiment on DNA con- 
fined between two parallel plates [IQ] , the average radius 
of gyration was measured from the configurations pro- 
jected onto a plane parallel to the confining glass plates. 
Ref. [l^l showed that with reducing plate-separation D, 
this projected radius of gyration first shrank, before ex- 
panding in accordance with de Gennes scaling £)~^/^. 
Similar non-monotonic features were also reported in 
Ref. This behavior should not be confused with the 
non-monotonicity in (-Rtot) (FigEl) and is clearly not in 
agreement with our assertion mentioned above. Ref. (ll| 
speculated that this non-monotonicity may be due to at- 
tractive interaction with the confining walls, as the de- 
tailed nature of polymer-wall interaction is both hard 
to control and to determine in experiments. It should 
be noted, however, that the non-monotonicity in the ra- 
dius of gyration calculated from the full 3d configura- 
tions of a confined self-avoiding polymer, as observed in 



our Langevin dynamics simulations (data not shown) and 
previous Monte-Carlo simulations [l^], has a different ori- 
gin. It is purely an entropic effect caused by the hard con- 
finement and, as we have shown in Sec. |lll is already to 
be expected on the basis of de Gennes scaling. Using our 
Langevin dynamics simulations, we have also evaluated 
the radius of gyration of the excluded volume polymer 
confined between two parallel plates from configurations 
projected onto a plane parallel to the plates, exactly in 
the same manner as described in Ref. These simu- 
lations showed a clear monotonic increase as D-^/^ with 
reducing D (data not shown}- This is in contrast to the 
observed behavior in Refs. [13, [HI, and in agreement with 
our assertion that polymer size in the unconfined direc- 
tions can only monotonically grow with increasing degree 
of confinement. 

Finally, the MET approach used in this paper can be 
easily extended to other confining geometries, e.g., cu- 
bic or cylindrical pores. The cylindrical confinement is a 
natural choice for the study of bacterial chromosomes [H) ■ 
Within a biological cell, binding of proteins on DNA may 
enhance (or reduce) the effective persistence length of 
DNA. Thus it would be interesting to examine the ef- 
fect of varying bending stiffness on confined polymers. 
A similar MET method may potentially also be used to 
describe tethered polymers. For this case one may need 
to self-consistently determine the step-size of the effec- 
tive Gaussian polymer as a function of distance from the 
tethered end. Thus, this mean field approach has the 
potential to provide analytical access to many situations 
of biological relevance. 
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Appendix A: Calculating Green's function 

A separation of variable Gq = /(L)?/'(r||)^(z) in 
Eq. (Uni) leads to, 

fdL 6 i/) " ^ 6 ^ ' 

where E is an arbitrary constant. Further, let us assume. 



dz^ 

2 , 1,2 



(z) = 



V|+k|U(r||) = 0. 



(A2) 
(A3) 
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The Dirichlet boundary condition Go{z = 
0,D) = requires the solution, (t)n{z) — 
y^2fDsm{k^z) ,/2/Dsm{k^z') with fcj = (riTr/D), 
n being a positive integer. The solution in the x and y 
directions gives = exp(«k|| .r||) with a continuum 

of k| | mo des in presence of an open boundary condition. 
Eq. (|A1I) leads to the identity, 




Thus the solution of the differential equation Eq. (|20)) has 
the form exp(— exp(ikj[.r||) sin(A:"z) sin(fc"2;'). The 
general solution is obtained by summing (integrating) 
over all possible modes n (k||). 



Go 



dkii 



D J (27r)2 

^ ' n=l 



-((|| /6)k^L-(ii/6)(W-D)'i *k|| .(r,, -r' 



II' sm y'jj'^ J sm | -^z' 



27r/||L 



'3('-||-'-l|)'/2i||i 



71 TT 



UTT 



(A5) 



Appendix B: Contributions from excluded volume 
interactions 



In this appendix we calculate S^^ and S^. 
The first order contributions due to the excluded vol- 
ume term in the Hamiltonian are of the form 



(A(r(0),r(L))/3A%)o = 

^1 
Z2 

r(L) 

V[r{s)]S (r(s) - r(s')) exp {-/3Ho[r(s)]} 



^-b J dr{0) J dr{L)A{r{0),r{L)) ds ds' 



(0) 

^^\b J dr{0) J dr{L)A{r{Q),r{L)) 

ds ds'£;(s,s'|r(0),r(L)), 
Jo Jq 

where the integration over end-points are performed for 
all possible r(0) and r(L). We consider the evaluation 
of the kernel E (s, s'|r (0) , r (L)) of the inner integration, 
assuming for the moment that s > s' . The Chapman- 
Kolmogorov property 

G(r,L|r',0) = J dr"G (r,i|r", L - s) G (r",s|r', 0) 

(Bl) 

then allows us to write 

E{s,s' I r(0),r(L)) = J dr'Go(r(L),L | r',s) 

Go(r',s I r',s')Go(r',s' | r(0),0). 

(B2) 

We also use the fact that since the Green's function fac- 
torizes as Gq = G||Gj_, the separation of variable works 



for E^E^i{s,s' I ^{0),^{L))xE^{s,s' \ r^{0),r^{L)) 
as well. 

As the longitudinal part is unaffected by the con- 
straints, it is the simplest case. Indeed, as 



G|| (r',s|r',s') 



27r?|| (s- s') 



(B3) 



does not depend on the intermediate position r'| we im- 
mediately find that 

£;|j(5,s' |r||(0),r||(L)) 

^ -G||(r||(L),L-(.-,s')|r||(0),0) 



277^11 (S - S') 



(B4) 



A similar simplification does not hold for the transverse 
component as 



(B5) 

(with K = {l±/6)(Tr/ D)"^) does depend on r'j_. We there- 
fore find that 



Ej_{s,s' \r^{Q),r^{L)) 
1 



1 X — ^ r 2 r f2 2\ '/ 2 f2\ 



TiL ,n' ,nQ — l 

X [26nQ,nL ^n',ino+nL)/2 — Sn',{no-nL)/2 — Sn',{nL-no)/2] 

X sin (^r^(i)) sin (^^±(0)) , (B6) 
Again, we can write 

{RpAnb)o = ^'^\^ "^''Z '^^'Z '^^(^^ / ^^^-^^ 
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[r||(i:)-r||(Or^i?(s,5'|r(0),r(L)) 
Z-'^h I ds I ds' 



Thus, 



"'0 



dr||(0) J d^{L) [r||(L)-ri|(0)]2 i;,, 
J dr^iO) J dr^{L)E^ 

Z-^\h j\s'll{s,s')ll{s,s'). 

(B7) 



Similarly, 

{Rl(3Anb)o = Z-^h j ds I ds'll{s,s')ll{s,s') 
^ Jo Jo 

(B8) 



{pAnb)o = Z-'^b j\s j\s'll{s,s')ll{s,s'). 

(B9) 



In the above, we used the definitions 



Il{s,s') = j dr|,(0) j dv\\{L) E\^{s,s' \ r,, (0), r,, (L)) 

(BIO) 

Il{s,s') = j dr||(0) j d^{L) I r||(0),ri|(i)) 

X (r||(L)-r||(0)f (Bll) 
Il{s,s') = J dr^iO) J dr^{L) E^{s,s' | r^{0),r^{L)) 

(B12) 

Il{s,s') = J dr^iO) J dr±{L) E^{s,s' \ r^{0),r^{L)) 
X K(L)-r^(0))2. (B13) 



Wc can now evaluate the integrals over the end points. 
The longitudinal ones are easier to compute. 



27rZ||(s- s') 

^i|(^'^') = ^'2;^^iK737)4'"(^^ 

(B14) 

The only r^(0) and r±(L) dependent term present 
in E±{s,s' I r^{0),r^{L)) is q{r^{L),r^iO)) = 
sm{nLTrr±{L)/D) x sin(no7rr_L(0)/£)). To evaluate the 
integration in Ij_{s, s') we use the identity 



ni^,n' ,no — 1 
X [2i5no,ni, + '^n',(no+"i,)/2 
,(no— nL)/2 "n',{nL — 

no)/2\Jo{nL,no) 

(B16) 

where 



= ^^^^^ ~ - (-in- (B17) 



To calculate l]_{s, s') one requires to use the integration 
of sine function with powers, such as, dr sm{mrr/D), 

drr sm.{mrr/D), and drr'^ sm{n'Kr/D). This 
gives us, 



nj:,n',no— 1 
X [2^no,nL + ^n',{no+nL)/2 
— ^n',{no-nL)/2 — ^n' ,{nL-no) /2] 

X J{nL,no), (B18) 

where. 



x [i-(-ir 



riLnoTr^ 



:{-ir'{-iy 



1 



^[(-ir.g.^+2(i-(-ir)] 



X [l-(-l)"'^]. 



(B19) 



We used s > s' in the above. Now, using the identity 

-I pL pL pL pL 

- ds ds' = ds' / ds (B20) 
2 Jo Jo Jo J s' 



we have 



5\ = Z-^b f ds' f ds Il{s,s')[ll{s,s') - {R\)oll{s,s')] 

Jo J s' 

(B21) 

5\ = Z-'bj ds' j ^ ds ll{s,s')[ll{s,s')-{Rl)oll{s,s')]. 

(B22) 

We now have all the ingredients to calculate the pertur- 
bation corrections. 
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1. Parallel Component S^, 



We have, 

|Jf(s,s')-(Ri)„7f(s,s')l 

j)x(|¥i-(»-.'))-|'„i) 



3W^ 



27r^|| (s — s' 



TT 



(B23) 

Thus m evaluating the integral in (Jy (Eg. IB2l|) . we need 
to do the following integration, 

M = e"*^"' [ ds' [ ds e-''^(""-"i)e-'^'*'("o-"") 

Jo Js' 

(B24) 

where the only function of (s, s') in the integrand conies 
from Ij_{s, s'). 
Thus, we obtain, 

6 1°° 
S\ = ~z^^- E M{nL,n',no;K) 

n L ,^0 — 1 
[2<^rio,ni + '^n%(no+ni)/2 

i.e., Eq.[36l 



2. Perpendicular Component 



The only function of (s, s') in Eq. lB22l is exp[— /ts(ri'^ — 
n\)] exp[— Ks'(ng — n'^)] due to the function [/^(s, s') — 
{R\)qI^^_{s,s')] and l/(s - s') due to I|^(s,s')- Thus we 
have to evaluate the integration, 



K 



^—K,Lr 



L pL 

ds' I ds 



/ '2 2 \ 1(2 I2\ 

^ — Ks{n —rij^j^ — KS [riQ—n ) 



{s - s') 



(B25) 

The term by term calculation of this quantity for each 
(tt-Oj n' ,ni) encounters a pole at s = s'. To avoid the pole, 
the integration to evaluate Eq. (jB25p can be rewritten as, 

2 

K(nL,n',no;K) = e""^"^ lim / ds' / ds 

2 ^,2 



Js' 

(11 2 \ ' / 2 /2 \ 



{s - s') + a 



(B26) 



Note that this is equivalent to replacing the delta- 
function overlap-interaction by its Gaussian representa- 
tion but keeping a non-zero variance . Thus we obtain 
the expression (Eq. [ 



I n/,,n',no — 1 



Sn' ,(no-nL)/2 — ^n' ,{nL-no) /2] 



J{nL,no) 



{Rl)o 

£|2 



Appendix C: Bulk limit of D oo 

In this appendix we calculate and Sj^ in the bulk 

limit of -D — ?> oo. 

It is clear that, since in this limit G± itself is indepen- 
dent of D, and 6^ are also independent of D. Using 
the continuum approach at Z? — >■ oo. 



2-kIj_{s — s 



G^{z{L),L~{s-s') I z(0),0) 



X exp 



2ttI^{s-s') y 2ttI^{L- {s- s')) 
3 (z(i)-z(0))2 1 



2lj_{L- (s- s')) 
Then using Eq.[B12]we get 



(CI) 



Ilis.s') = J dziO) J dz{L)E^is,s' I r^{0),r^{L)) 

Remember that, now, the orthogonal component of par- 
tition function has also got redefined to Zj_ = D. Then 
using Eq. IBITI 

= Z-^h ds' ds ll{s,s')[ll{s,s')-{R])Ql1{s,i 
= {W^D)-^h{~W^/TT) 

X dJ^— f ds' [ ds 

b ~ 



{s^s') 



X 1x3/2. 

2ttI±_ 3 



(C 



where we also used Egs. IB14l and[?fl Using Eg. IB13I we 
find 



Il{s,s') 



dz{0) / dz{L). 



2ttIj_{s - s') 



2ttI±{L^{s - s')) 
' LiL-is-s')) 



y 2ttIj_{s- s') 
and using Eq. [51] we get 



D- 



(C4) 



r dzT dz'iz-z'fe-'-''''^^ 

J— oo J— oo V -' 



f°° dz f°° dz't 
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h_L 



(C5) 



Therefore, 



Il{s,s') - {RDollis^s') = -D^{s-s'f'\ (C6) 



and 



I^,{s,s')[ll{s,s')-{Rl),ll{s,s')] 
3 „ 17 



27r/||(s - s') V 67r 



/Nl/2 



1/2 



-{s-s')-"'. 



(C7) 



Since in this hmit Z = W^D, using Eg. IB22l we find, 
5\ = Z-^h j ds' J ds I^\{s,s') 

2 ('„ „'^ / d2 \ tO f „ „'\ 



X [liis,s')-{Ri)olUs,s')] 



hi 



1/2 



.^3/2 



(C8) 



Similarly, using Eq. IC5I in Eq. [SU we obtain 



(C9) 



where Zod = . Thus, 



A'' - 

"2ci — 



where ho plays the role of interaction strength. Thus, 
^11 = 5}, +5\^Q implies 



Ll^d ^ ^0 j-2 

; 27r 

or, - m'/'L'/' 



which leads to the scaling form. 



(Dl) 



We expect that the polymer would behave like a 
pure 2d polymer in the limit of extremely small plate- 
separation D. This limit is achieved when = {R\^^- 
We have seen that the strength of the interaction h in 3d 
has the dimension of length, whereas the same strength 60 
is dimension-less in 2d. This shows that the only intrinsic 
microscopic length scale in the system is segment-length 
I. Thus we need to express h — b^l to search for the pure 
2d limit of confined systems. With this substitution, 



Appendix D: 2d limit of confined system 

In pure 2d, the perturbative contribution due to the 
inter-segment interactions takes the form, 

^Id = Z^dbo f ds' f ds[ll{s, s') ~ (i?|)o/|?(s, 3')] 

Jo J s' 



(i?f) ^ {holf'L^'^ 



D 



1/2 



(D2) 



Therefore, equates (i?2d) when D ^ I, i.e., when 

the plate-separation D becomes as small as the polymer 
segment-length I. 
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